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ABSTRACT 

Echelle spectroscopy of H 2 2.122 p,m, [Fell] 1.644p,m and Bry line emission from a 
very young planetary nebula, IRAS 21282+5050, and from four proto-planetary neb¬ 
ulae, IRAS 19343+2926 (Ml-92), IRAS 17150-3224 (AFGL 6815), IRAS 17423-1755 
(Hen 3-1475), and IRAS 17441-2411 is presented. H 2 line emission is detected in dis¬ 
crete shock fronts in the lobes of each nebula, regardless of source spectral type (al¬ 
though non-detections in IRAS 09371+1212 (Frosty Leo) support claims that late 
spectral types do not produce bright H 2 line emission). In IRAS 17150-3224 we also 
uncover possible signs of rotation, as would be expected if the H 2 features were ex¬ 
cited in a magneto-centrifugal disk wind. [Fell] 1.644^m emission was detected in 
only one source, Ml-92 (notably, the source with the brightest H 2 features). Again, 
the emission is predominantly excited in high-velocity shocks in the bipolar lobes of 
the PPN. The H 2 and [Fell] observations of M 1-92, and the complex H 2 profiles in 
IRAS 21282+5050, are explained using the shock models of Smith and collaborators. 
We show that bow shock models are generally able to account for the observed line 
profiles, peak velocities, the double-peaked profiles in IRAS 21282+5050 and the spa¬ 
tial distribution of H 2 and [Fell] in M 1-92. J-type bow models are adequate in each 
case, i.e. a strongly-magnetised wind is not required. 

Finally, Bry is detected in each of the five targets; in absorption in the G-type 
PPN though in emission in the O and B-type sources. Bry emission is detected pre¬ 
dominantly toward the near-lR continuum peak in each PPN, with only very weak 
emission detected in the extended lobes of the O- and B-type sources. In Bry, low peak 
velocities, though very broad profile widths, are measured in each target, regardless 
of nebula inclination angle. The emission must therefore derive from ionised regions 
in a fast wind very close to the central star (rather than from shocks in the bipolar 
lobes), or, in the late-type sources, from absorption in an equatorial torus. 

Key words: circumstellar matter - infrared: stars - infrared: planetary nebulae - 
ISM: jets and outflows - ISM: kinematics and dynamics 


1 INTRODUCTION 

The proto-planetary nebula (PPN) phase is a short-lived 
episode in the rapid evolution of stars from the asymptotic 
giant branch (AGB) to the planetary nebula (PN) phase 
(Kwok 1993). The central stars of PPN are still too cool 
to ionise the slow-moving circumstellar shell ejected during 
the AGB phase, yet they appear to drive collimated, high- 
velocity winds which shock and shape this shell, and almost 
certainly entrain slow-moving AGB ejecta to produce a fast 
molecular wind (Lee & Sahai 2003). 

From recent observations, particularly high-resolution 


imaging studies (e.g. Sahai et al. 1998; Bujarrabal et al. 
1998a; Sahai, Bujarrabal & Zilstra 1999; Ueta, Meixner & 
Bobrowski 2000; Su, Hrivnak & Kwok 2001), it has become 
clear that the rapidly-evolving PPN phase ultimately shapes 
the morphology of the resulting PN. At some point during or 
soon after the detachment of the AGB envelope, the mass- 
loss geometry changes from more-or-less spherically sym¬ 
metric to axially symmetric. Bipolar morphologies, knotty 
jets and Herbig-Haro (HH)-like “bow shocks” are all mani¬ 
festations of this evolutionary process. Indeed, they seem to 
be present even in relatively “young” PPN. 

Recently, Bujarrabal et al. (2001) have found that 
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radiatively-driven winds cannot account for the high mo¬ 
mentum and energy implied by their CO observations of 
PPN fast winds. This has prompted theorists to look be¬ 
yond the classical “interacting stellar winds” model of Kwok 
et al. (1978) and investigate whether an accretion disk sce¬ 
nario, similar to that used to explain jets from young stars 
and active galactic nuclei, could be used to explain both the 
high degree of collimation and the point symmetry seen in 
many PPN jets (Soker & Rappaport 2000; Frank & Black- 
mann 2004). In these models, the accretion disk must form 
through binary interactions. Magneto-centrifugal launching 
from the surface of the disk is then a means of convert¬ 
ing gravitational energy into kinetic energy in a fast wind. 
If magneto-hydrodynamics (MHD) do dictate the energet¬ 
ics and collimation of PPN outflows (see also Garci'a-Segura 
et al. 1999; Blackman et al. 2001), then they will also de¬ 
fine the physics of the shocks in these flows, an effect which 
can be investigated with high-resolution observations of the 
emission regions associated with the shocks. 

We therefore observed a modest number of PPN and 
young PN at high spectral resolution. Our sample in¬ 
cludes reddened PPN (IRAS 17150-3224, IRAS 17441-2411, 
IRAS 09371-bl212 [Frosty Leo] and IRAS 17106-3046) 
as well as more evolved PPN and young PN (M 1-92, 
IRAS 17423-1755 and IRAS21282-1-5050), sources with ear¬ 
lier spectral types that are on the brink of the transition 
to the PN stage. We sought observations in a number of 
atomic and molecular gas tracers. H 2 l-OS(l) emission at 
2.122 pm is an important probe of excited molecular gas in 
the detached AGB envelope (>90% of the AGB ejecta will be 
atomic/molecular), while Bry recombination line emission is 
potentially a useful tracer of the fast, ionised wind. [Fell] at 
1.64 pm could also be an excellent tracer of moderate-to-high 
excitation shocks in the fast wind and/or the interaction 
zone between the wind and AGB envelope (the Fed- ioni¬ 
sation fraction peaks at T ^ 14,000 K for uh ~ 10®cm“®; 
[Fell] has been shown to be a useful diagnostic in studies 
of HH shocks, e.g. Nisini et al. 2002; Davis et al. 2003a). 
In the dense, inner-core regions of each PPN, these near-IR 
lines are also less affected by extinction than their optical 
counterparts, and so are expected to be ideal probes of the 
dynamical interaction between the fast wind and AGB en¬ 
velope. 


2 OBSERVATIONS AND DATA REDUCTION 

Echelle spectra at 1.6440 pm, 2.1218 pm and 2.1662 pm were 
obtained at the 3.8 m United Kingdom Infrared Telescope 
(UKIRT), on Manna Kea, Hawaii, on various dates in 
Spring 2004 (see Table 1). The cooled grating spectrome¬ 
ter GGS 4 (Mountain et al. 1990) was used. This instru¬ 
ment is equipped with a 256x256 pixel InSb array and a 
31 lines/mm echelle grating. The pixel scale is 0.4l"x0.88” 
(0.41" in the dispersion direction). A 2-pixel-wide slit was 
used, resulting in a velocity resolution of ~ 16 km s“^ (this 
being the Full Width at Half Maximum (FWHM) of unre¬ 
solved sky or arc lamp lines). 

At the above wavelengths the H 2 l-OS(l) (Avac = 
2.1218334 pm; Bragg, Smith & Brault 1982), Hi Bry 
(Avac = 2.166167 pm) and [Fell] DDt /2 — a‘^Fg /2 (Avac = 
1.643998 pm; Johansson 1978) transitions are well centred 


on the array. The wavelength coverage at each of these wave¬ 
lengths was approximately 2.115-2.128 pm, 2.158-2.173 pm 
and 1.638-1.649 pm (these ranges shifted slightly from night- 
to-night, depending on the setting of the echelle grating an¬ 
gle). 

For each target, and for each line, an exposure at a 
“blank sky” position was followed by hve on-source “ob¬ 
ject” exposures, the former being subtracted from all of the 
object frames. This sequence of six frames was repeated two 
or three times in each line to build up signal-to-noise. To 
avoid saturation on the bright continuum associated with 
each PPN, relatively short exposures were used; 60sec (x2 
coadds) for the two K-band lines, and 80 sec (x 2 coadds) in 
the H-band. The total integration time on each source was: 
30 mins in H 2 - central slit (in all seven targets); 20 mins 
in H 2 - left and right slits (in the five targets with H 2 de¬ 
tections); 26.5 mins in [Fell] (40 mins for IRAS 17150-3224 
and IRAS 17423-1755), and 30 mins in Bry. 

With the above exposure times, at the high spectral res¬ 
olution of the instrument {R ~ 18, 750), there is essentially 
no detectable thermal emission (or associated shot noise) 
from the sky or background between the well-resolved OH 
sky lines. Thus, the process of “sky subtraction” mainly 
serves to remove the OH lines and any “warm” pixels that 
remain after bad-pixel masking and flat-helding (all data 
frames are bad-pixel-masked, using a mask taken at the start 
of each night, and flat-fielded, using an observation of an 
internal blackbody lamp obtained before each set of target 
observations). At 2.122 pm and 2.166 pm the H 2 l-OS(l) and 
Bry lines are well separated from any bright OH sky lines. 
Thus, any residual sky lines in the reduced spectral images 
resulting from imperfect sky subtraction do not affect the 
data. However, the [Fell] 1.644 pm line is close to an OH 
line, the 1.644216 pm 5-3 R2 transition, so residual emis¬ 
sion from this OH line was removed using the STARLINK 
package FIGARO, specifically the polysky routine. This rou¬ 
tine fits a polynomial to regions above and below the target 
spectrum (along columns in the spatial/slit direction): the 
ht to these sky regions is then used to subtract sky emission, 
column-by-column, across the spectral image. 

Goadded H 2 and [Fell] spectral images were wavelength 
calibrated using OH sky lines, there being four lines well 
spaced across the array in both cases (Rousselot et al. 2000). 
The hrst raw object frame observed at each wavelength on 
each source was used as the reference frame. The STAR- 
LINK.FIGARO software packages used to do this also cor¬ 
rect for distortion along the columns in each image (i.e. along 
arc or sky lines), via a polynomial ht to the OH lines in each 
row. Examination of the OH lines in distortion-corrected and 
wavelength-calibrated raw frames showed that the velocity 
calibration along the length of the slit, measured from Gaus¬ 
sian hts to the sky lines, is extremely good for the H 2 and 
[Fell] data. The relative velocity calibration along the slit 
and between adjacent slits and targets is estimated to be 
accurate to ~ 3km s“^ . 

Unfortunately, only one bright OH line was evident 
in the raw Bry spectral images, at 2.1711170 pm. A pair 
of weaker lines appears in some the data at the blue end 
of the spectrum, and a second pair of weak lines occurs 
near the centre of the array. However, these lines were gen¬ 
erally too faint for accurate calibration, and were in any 
case unresolved by Rousselot et al. (2000). We therefore 
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Table 1. UKIRT/CGS4 echelle spectroscopy observing log 


IRAS 

number 

Other 

name 

Central^ 

star 

(kpc) 

RA® 

(2000.0) 

Dec® 

(2000.0) 

Date'^ 

(H 2 ) 

Date^ 

(Br7) 

Date'^ 

([Fell]) 

Slit® 

P.A. 

e® 

ReP 

21282+5050 


09 

3 

21 29 58.5 

51 04 00 

040728 

040728 

040728 

160° 

0° 

1 

19343+2926 

M 1-92 

B0.5V 

2.5 

19 36 18.9 

29 32 51 

040430,040705 

040530 

040530 

133° 

35° 

2 

17423-1755 

Hen 3-1475 

Be 

5 

17 45 14.2 

-17 56 47 

040615,040705 

040624 

040624 

135° 

0 

0 

10 

3 

17441-2411 

AFGL 5385 

GO 

4 

17 47 13.5 

-24 12 51 

040625,040703 

040625 

040705 

197° 

<30° 

4 

17150-3224 

AFGL 6815 

G2I 

2.4 

17 18 19.8 

-32 27 21 

040615,040705 

040625 

040625 

125° 

<30° 

5 

09371+1212 

Frosty Leo 

K7H 

3 

09 39 53.6 

11 58 54 

040428 

- 

- 

156° 

15° 

6 

17106-3046 


? 

7 

17 13 51.6 

-30 49 40 

040625 

- 

- 

128° 

52° 

7 


^ Optical spectral type of the central object 
^ Estimated distance to the source 
® Central source position. 

UT date of the H 2 2.122/im, Bry and [Fell] 1.644/rm observations of each target (year/month/day). The “left” and “right” H 2 slits 
referred to in Section 3 were, in most cases, observed on the second date in column 7. 

® Position angle, measured East of North, of the echelle slit. The same angle was used at each wavelength, and for each of the three 
adjacent H 2 slit positions (the angle corresponds to the axis of the bipolar lobes in each PPN). 

® Inclination angle of the bipolar nebula with respect to the plane of the sky. 

^ References to D and 0: 1) Meixner et al. 1997; 2) Bujarrabal et al. 1997, Bujarrabal et al. 1998b; 3) Riera et al. 1995, Borkowski & 
Harrington 2001; 4) Su et al. 1998, Weintraub et al. 1998; 5) Loup et al. 1993; 6) Mauron et al. 1989, Roddier et al. 1995; 7) Kwok et 
al. 2000 


used telluric absorption features in the standard star spec¬ 
trum to wavelength calibrate the central region in each 
Br 7 spectral image. Note that in this case no correction 
was made for optical distortions along the slit axis. Fortu¬ 
nately, the Bry emission was confined to a region very close 
to the source continuum in each spectral image, so distor¬ 
tion was not considered to be a major source of uncertainty 
in the wavelength/velocity calibration. Three or four tel¬ 
luric lines were observed in each standard star spectrum; 
these were identihed in the Arcturus atlas of Hinkle, Wal¬ 
lace & Livingstone (1995) and their wavelengths measured 
from the solar atlas available from the BAse de donnees 
Solaire Sol 2000 solar archive hosted by I’Observatoire de 
Paris (http://bass2000.obspm.fr/solar_spect.phpl. Bry ab¬ 
sorption lines, seen in each standard, were also used (after 
correction for the radial velocity of the standard). The OH 
lines noted above were subsequently used to check the cal¬ 
ibration. The telluric features were relatively weak in most 
cases, and their wavelengths only measured to an accuracy 
of ~ 0.1 Angstrom. The velocity calibration in the Bry data 
is therefore only considered to be accurate to ~ 5 — 10 km 
s-i. 

At all three wavelengths velocities were calibrated with 
respect to the kinematic Local Standard of Rest [LSR]. They 
were not corrected for the systemic velocity in each region. 

Finally, spectra of A or G-type standard stars were ob¬ 
tained and reduced in a similar fashion. In each case an ex¬ 
tracted, wavelength-calibrated spectrum was “grown” into a 
2-D image so that division by this image would correct each 
reduced PPN spectral image for telluric absorption (Bry ab¬ 
sorption lines were removed from the standard star data 
before they were used to calibrate the 2.166 /rm images). Al¬ 
though most of the data were obtained in photometric con¬ 
ditions, variable cirrus on 15 June, 24 June and 3 July intro¬ 
duced 0.5-1.5 magnitudes of attenuation. The absolute flux 
calibration of these observations (i.e. the centre H 2 slit in 



Figure 1. Narrow-band 2.122/rm image of IRAS 21282+5050, 
acquired to establish the orientation of the bipolar lobes in this 
young PN. Note that the image is not continuum-subtracted; nor 
is it flux calibrated, so the contours are arbitrary (although they 
approximately represent 1,2,4,8,16,32,64,96,128,192,320 and 480- 
times the standard deviation in the blank-sky around the source). 
North is up and east is to the left. 

IRAS 17150-3224 and IRAS 17423-1755, the [Fell] and Bry 
observations of IRAS 17423-1755, and the left and right H 2 
slits in IRAS 17441-2411) will therefore only be accurate to 
within a factor of 2-3. The seeing was in the range 0.5”-0.8”, 
i.e. equal to or somewhat less than the slit width and offsets 
used for the adjacent H 2 observations. 


3 RESULTS 

In this paper we have targeted a very young PN and six 
PPN; the sources are listed in order of spectral type in 
Table 1. The spectral energy distribution (SED) of each 
source is double-peaked, though in most cases dominated 
by a strong infrared component, as expected for a detached. 
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IRAS21 282+5050(cen) IRAS21282+5050(cen) 





Figure 2. PV diagrams showing the H 2 2.122/im emission in IRAS 21282+5050. The top of each spectral image is orientated roughly 
to the north (p.a. ~ —20°). With respect to the centre slit (cen), the “left” and “right” spectra were offset 1.2" to the east and west, 
respectively. The data from the central slit are shown before and after continuum-subtraction (top two panels); only continuum-subtracted 
PV diagrams are shown for the left and right slits (bottom two panels). The contour levels measure 0.005,0.015,0.03,0.06,0.09,0.12,0.24,0.48 
Jy per arcsec along the slit in all four plots. The radial velocities of emission peaks, measured from Gaussian fits to spectra extracted 
over 3 rows (2.7") centred on the peak, are labeled. The velocity zero point is not corrrected for the systemic LSR velocity. 


dusty shell surrounding a stellar photosphere (e.g. Kwok 
et al. 1996; Ueta et al. 2000). Bipolar lobes are evident in 
published, high-resolution images of each target (except for 
IRAS 21282+5050, though see Figure 1). The lobes in each 
case extend over at least a few arcseconds in projection on 
the sky. In the optical, the aspect ratio of the bipolar lobes 
(length/width) is typically 2-5. 

For each source the spectrometer slit was centred on 
the bright near-IR continuum peak; the slit was aligned 
along the axis through the bipolar lobes, as defined by 
high-resolution images (e.g. M1-92 - Bujarrabal et al. 
1998a; IRAS 17423-1755 - Reira et al. 2003; IRAS 17441- 
2411 - Su et al. 1998; IRAS 17150-3224 - Kwok et al. 1998; 
IRAS 09371+1212 - Sahai et al. 2000; IRAS 17106-3046 - 
Kwok, Hrivnak & Su 2000). We refer to this position as the 
“central slit” (cen) position in each case. We found no high- 
resolution narrow-band images of IRAS 21282+5050 in the 
literature (although see the K-band image of Latter et al. 
1995). We therefore obtained a narrow-band (A/dA ^ 100) 
2.122 fim image of this target under good seeing conditions 
(^ 0.6”) with the UKIRT Fast Track Imager (UFTI) (Roche 
et al. 2003). An adjacent continuum image was not obtained; 
nor were the data flux-calibrated. In this image (Figure 
1), IRAS 21282+5050 was seen to be slightly extended in 
a north-north-west to south-south-east direction, as seems 
to be the case in the K-band image of Latter et al. (1995): 
a position angle (p.a.) of 160° was thus adopted for this 
target. 

Although H 2 observations were obtained for seven tar¬ 


gets, line emission was detected in only the hrst five sources 
in Table 1. The upper limit to the line flux in the lobes 
of the two non-detections was ^ 5 mjy arcsec”^ (~ 1 x 
10”^® W m~^ arcsec”^). When H 2 was detected in the cen¬ 
tral slit position, i.e. along the polar axis of a source, H 2 
data at adjacent, parallel slit positions were also obtained. 
In each case the slit was offset ^ 1" either side of the central 
axis. The exact offset - reported in the H 2 position-velocity 
(PV) diagram Figure captions - is measured perpendicular 
to the slit axis. We refer to these offset positions as “left” and 
“right”. In this way, H 2 spectra were obtained at three ad¬ 
jacent, parallel slits across the width of the hrst hve sources 
listed in Table 1. 

[Fell] 1.644 fim and Bry observations were also then ob¬ 
tained for these hve targets, with the slit centred on the 
same continuum peak position and using the same posi¬ 
tion angle. In [Fell] and Bry only this one central slit po¬ 
sition was observed in each source. Bry was detected in 
all hve targets. However, [Fell] was only detected in M 1- 
92. The upper limit to the [Fell] hux in the lobes of the 
IRAS 17150-3224, IRAS 17423-1755, IRAS 17441-2411 and 
IRAS 21282+5050 nebulae was ~ 3 — 10 mJy arcsec“^ 
(^ 0.5 — 2.0 X lO”*^® W m“^ arcsec"^). 

Below we discuss the sources individually. H 2 PV dia¬ 
grams are shown for each target. For the central slit position 
a PV diagram showing the line and continuum emission is 
presented along-side a continuum-subtracted PV diagram. 
In the latter the continuum has been removed from the spec¬ 
tral image, row-by-row, by htting the continuum on either 
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side of the line emission region with a second-order polyno¬ 
mial. The PV diagram showing both the line emission and 
the continuum serves to illustrate the relative strength of 
the two emission components. For the left and right slit po¬ 
sitions, only continuum-subtracted PV plots are presented. 
In all PV diagrams positive offsets are approximately to the 
north or north-west, depending on the slit position angle. 
Slit angles are listed in Table 1 and in the Figure captions. 

3.1 IRAS 21282+5050 

IRAS 21282+5050 is probably slightly more evolved than 
the other targets discussed in this paper. It is considered 
to be a very young PN with a WR central star (reclassified 
as 09 by Crowther, De Marco & Barlow 1998) surrounded 
by a dusty core with a central hole (Meixner et al. 1997; 
1998). It exhibits weak radio continuum emission at 2 cm 
and 6 cm (Likkel et al. 1994) and has probably evolved from 
a carbon-rich, intermediate mass star. 

Our narrow-band, 2.12 image of IRAS 21282+5050 
(Figure 1) shows the target to be elongated in two direc¬ 
tions. The outer, faint contours extend toward the north- 
north-west and south-south-east, along an axis orientated 
at a position angle of about 160°. We assume this to be the 
orientation of the bipolar lobes and adopt this angle for our 
spectroscopy. However, the brighter contours in the core of 
IRAS 21282+5050 are orientated perpendicular to this axis, 
at a position angle of 75°(±10°). The near-IR image in Fig¬ 
ure 1 bears a striking resemblance to the 12.5 ^m image 
of Meixner et al. (1997). They observe the same elonga¬ 
tion toward the north-north-west and south-south-east (p.a. 
~165°) and also, in their deconvolved image, resolve the in¬ 
ner core into two peaks, aligned along an axis orthogonal 
to the fainter emission (p.a. ~70°). They interpret their im¬ 
ages and the source’s spectral energy distribution in terms 
of dust emission from a disk viewed edge-on, implying that 
the bipolar lobes must lie close to the plane of the sky. 

Hora et al. (1999) included IRAS 21282+5050 in their 
extensive spectroscopic survey of PN. They obtained spectra 
at two positions, toward the core and toward a position offset 
3”N, 3”W. They detected H 2 emission lines from the latter. 
Davis et al (2003b) later presented an integrated spectrum 
of this target; they too detected line emission from H 2 . 

In the central slit data in Figure |^we see H 2 2.122 /rm 
emission along the slit axis and toward the source contin¬ 
uum position. As we demonstrate in subsequent sections, 
this is somewhat unexpected, particularly for a source with 
bipolar lobes orientated close to the plane of the sky: in the 
other targets discussed later, the H 2 is confined to discrete 
knots in the bipolar lobes. The H 2 near the continuum po¬ 
sition in IRAS 21282+5050 is therefore unusual. It may be 
excited either very close to the central star, perhaps in the 
extended wings of bow shocks in the bipolar lobes, or in a 
second outflow component orientated along the line of sight 
and therefore in the equatorial plane. Equatorial H 2 fea¬ 
tures have been observed in at least one other PPN, namely 
CRL 2688, where multiple bow-shaped shock structures are 
spatially resolved in both the bipolar lobes and in the equa¬ 
torial regions (Sahai et al. 1998). 

Our three near-IR slits were orientated perpendicular to 
the dust ring found by Meixner et al. (1997). Consequently, 
the spatial variations along each slit probably correspond to 


IRAS21282 + 5050 (left/centre/right) 



LSR Velocity (km/s) 

Figure 3. H 2 spectra across the width of IRAS 21282-1-5050. The 
spectra represent three rows (~2.7”), centred on the source con¬ 
tinuum, extracted from the “left”, “centre” and “right” slit po¬ 
sitions (plotted with thick-dash, full, and thin-dash lines, respec¬ 
tively) . 

offsets along the bipolar lobes, and the emission above and 
below the continuum in each PV diagram must be associated 
with the bipolar lobes. In other words, each PV plot could 
include H 2 line emission from an equatorial flow and from 
the bipolar lobes. 

In the continuum-subtracted H 2 PV diagram for our 
central slit, the emission describes a ring in position-velocity 
space. From a slice taken along the slit axis (i.e. along a sin¬ 
gle 0.52 A/7.4 km s“^-wide column in the PV diagram), 
the diameter of this ring is found to be of the order of 
4.4”(±0.9”). This diameter probably represents the distance 
between emission features in the bipolar lobes, above and be¬ 
low the source continuum position. (It does not correspond 
to the diameter of the dust torus observed by Meixner et al. 
1997.) 

Continuum-subtracted PV diagrams from the “left” 
and “right” slits are also shown in Figure Spectra ex¬ 
tracted from the middle three rows (offset = 0”) in each 
of the spectral images are shown in Figure 0 The broad¬ 
est profile is seen in the centre slit data. Here the emission 
line is double peaked, and can be well fitted with a two- 
component Gaussian profile with peaks separated by 22 km 
s“^. The spectra from the left and right slit positions, al¬ 
though slightly skewed to blue and red-shifted velocities, 
peak at an intermediate velocity of ~ 16 km s“^. These spec¬ 
tra - all obtained toward the central continuum position - 
may therefore be understood in terms of: (1) emission ex¬ 
cited in the wings of extended bow shocks seen side-on (the 
double-peaked profiles representing emission from the near 
and far sides of the bow), or (2) an expanding shell or curved 
shock directed toward the observer. These two scenarios will 
be investigated in Sect. 4.3.1, where we discuss whether an 
equatorial wind is needed in addition to an orthogonal, bipo¬ 
lar flow to explain the H 2 profiles in Figures l^andO 

3.2 Ml-92 

M 1-92 (IRAS 19343+2926) is a bipolar, oxygen-rich PPN 
(spectral type B0.5V). The relatively-high temperature of 
Its central star (~ 20, 000 K) suggests that it is probably 
more evolved than the other PPN targets in Table 1. Bu- 
jarrabal, Alcolea & Neri (1998b) estimate a kinematic age of 
about 900 yrs for M 1-92. The nebula comprises two bipolar 
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Figure 4. PV diagrams showing the H 2 2.122 /jm emission in Ml-92 (see Figure 2 for details). The top of each spectral image is orientated 
roughly to the north-west (p.a. ~ —47°); the “left” and “right” spectra were offset 1.2” to the north-east and south-west, respectively. 
The contour levels measure 0.04,0.08,0.16,0.32,0.64,1.28,2.56 Jy per arcsec along the slit in all four plots. 


lobes of scattered, polarised light (Herbig 1975; Trammell 
& Goodrich 1996; Bujarrabal et al. 1998a). Spectra from 
these lobes consist of a stellar continuum, scattered permit¬ 
ted Balmer and Hel lines, and a variety of low-excitation for¬ 
bidden emission lines, the latter being conhned to fast (200- 
300 km s“^) shocks in the bipolar lobes (Solf 1994). The two 
lobes are separated by a dark lane, although the central star 
is just visible because of the inclination of the north-western, 
blue-shifted lobe toward the observer. In optical HST im¬ 
ages the continuum in the lobes appears smooth, while the 
shock-excited [SII] and [OI] is shown to delineate knotty, col¬ 
limated jets that are roughly, though not precisely, aligned 
with the axis of the nebula (Trammell & Goodrich 1996; 
Bujarrabal et al. 1998a). 

Two hollow molecular shells surround the bipolar lobes 
(Bujarrabal et al. 1997). Observed in GO, these probably 
represent the ejecta from the AGB phase; gas swept up by 
the fast, post-AGB wind. They contain most of the nebula 
material and are expanding at high speed (up to 70 km s“^). 
Bujarrabal et al. also report the detection of fast molecular 
clumps at the tips of the bipolar lobes/shells (i.e. along the 
jet axes), as well as a dense, expanding, equatorial torus, 
with a diameter of about 2”-3”. 

In near-IR images of M 1-92, knots of H 2 2.122 /rm emis¬ 
sion are observed along the axes of the two bipolar lobes (Bu¬ 
jarrabal et al. 1998a; Davis et al. 2003b). The emission from 
these regions is evident in our echelle data; spectroscopy at 
the three parallel slit positions across the width of M 1-92 
are shown in Figure^] Discrete H 2 emission features are seen 
in all three slits; H 2 is not detected toward the continuum 
position. The emission knot to the north-west is blue-shifted 
(radial velocity, Flsr ~ —30 km s“^) while the emission in 


the counter-lobe is red-shifted (Tlsr ~ +30 km s“^). When 
corrected for inclination, these values convert to a velocity 
along the bipolar axis of the system of about 50 km s“^, 
a speed that is comparable to the velocity of the molecular 
AGB ejecta (Bujarrabal et al. 1998b). 

Of the five sources detected in H 2 2.122 ^m emission, 
Ml-92 was the only target where [Fell] 1.644/rm was de¬ 
tected. In these [Fell] data (Figure we see blue-shifted 
and red-shifted features in the bipolar lobes similar to those 
seen in H 2 . Extracted H 2 and [Fell] spectra are compared in 
Figure]^ These data were extracted over the same regions 
in the two bipolar lobes. The [Fell] spectra peak at slightly 
more blue-shifted/red-shifted velocities (Flsr ~ —56 km 8“*^ 
and ~ 41km s“^); the [Fell] profiles are also considerably 
broader. Gaussian fits to the extracted H 2 profiles yield 
FWHM velocities of ^ 35 - 45 km s“^; the [Fell] spectra 
are roughly twice as wide, and they exhibit strong blue- and 
red-shifted line wings on the blue- and red-shifted peaks, re¬ 
spectively. At the 10% peak intensity level, the [Fell] lines 
are ~ 330 - 340 km s“^ wide. 

The blue and red [Fell] peaks are also slightly further 
apart in separation along the slit, as compared to their H 2 
counterparts. From profiles taken along the slit axis from 
the continuum-subtracted spectral images, integrated over 
the same velocity range in each case, the separation between 
the [Fell] peaks was found to be 5.53” compared to 5.20” for 
the H 2 peaks. An error of ^ 0.04” is assumed, based on the 
accuracy with which we were able to calibrate the pixel scale 
(from observations of standard stars); the errors from the 
gaussian fits to the peaks in the profiles were considerably 
smaller than this. 

Note, finally, that although no H 2 emission was detected 
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Figure 5. PV diagrams showing the [Fell] 1.644/rm emission 
through the centre of Ml-92. In the lower plot the continuum 
emission has been fitted and removed. The LSR velocity scale has 
been calibrated with respect to the wavelength of the [Fell] line. 
The contour levels measure 0.005,0.01,0.02,0.04,0.08,0.16,0.32 Jy 
per arcsec along the slit in both plots. 


Ml-92 



LSR Velocity (km/s) 

Figure 6. Plots comparing H 2 and [Fell] line profiles above and 
below the stellar continuum in Figure O averaged over 2.7" re¬ 
gions along the slit in both cases. The H 2 data are plotted with 
dashed lines, the [Fell] with solid lines. The H 2 spectra, extracted 
from the centre slit, have also been multiplied by a factor of 0.1. 


toward the central continuum position (0” offset in each of 
the PV diagrams), a spatially-compact feature is detected 
in [Fell]. In an extracted spectrum (not shown) this feature 
peaks at ~ —2 km s“^, i.e. roughly the nominal systemic 
velocity based on the blue- and red-shifted velocities seen 
in the lobes in H 2 and [Fell]. At only ~ 24km s“^ wide 
FWHM, this feature is also considerably narrower than the 
[Fell] profiles in the lobes, being only moderately broader 
than the instrumental profile. 


3.3 IRAS 17423-1755 

Optical images of the bipolar PPN IRAS 17423-1755 (Hen3- 
1475) reveal a spectacular S-shaped distribution of point- 
symmetric knots along the axis of the nebula, which extend 
over ~17" (Borkowski, Blondin & Hartigan 1997; Riera et 
al. 2003). Like M 1-92, this chain of knots is enveloped in a 
diffuse, dumbell-shaped, bipolar nebula. From optical spec¬ 
troscopy, velocities for the fast wind in the lobes of the neb¬ 
ula of ~ 1000 km s“^ have been inferred, and shocks are once 
again thought to play a major role in the excitation of the 
optical line-emission features (e.g. Riera et al. 2003). Riera 
et al. suggest that IRAS 17423-1755 may have evolved from 
a relatively high mass progenitor. VLA observations do re¬ 
veal a compact radio-continuum source near the central star 
(Knapp et al. 1995), so the central object may already be 
hot enough to produce an ionised region. 

H 2 observations at three positions across the width of 
IRAS 17423-1755 are shown in Figure|7| The outer slits co¬ 
incide with the edges of the diffuse lobes of scattered light, 
while the central slit covers many (though not all) of the op¬ 
tical knots along the central axis. Even so, H 2 line emission 
was only detected in the north-western outflow lobe (positive 
offsets in Figure O, probably toward the bright optical knot 
NWla (Riera et al. 2003). Our continuum-subtracted spec¬ 
tral images do not reveal emission coincident with the cen¬ 
tral star, although the shot noise associated with the bright 
continuum is much stronger than in the bipolar lobes. 

The H 2 emission in the north-western lobe peaks at an 
offset of ~3”; it is also strongest in our “right” slit position, 
which corresponds to the south-western edge of this lobe. 
The emission feature peaks at a highly blue-shifted veloc¬ 
ity (Vlsr ~ —70km s~ ). This is not too surprising when 
one considers that IRAS 17423-1755 is probably the most 
obliquely-orientated source in our target list (Table 1). The 
line profile is narrow (FWHM ~ 18km s“^) though there 
is also possibly a weak wing extended to higher blue-shifted 
velocities. 


3.4 IRAS 17441-2411 

With respect to angular size, aspect ratio and overall struc¬ 
ture, IRAS 17441-2411 (AFGL 5385) is very similar to 
IRAS 17150-3224 (discussed below). Both targets are also 
associated with late-type central sources. In HST images 
IRAS 17441-2411 comprises two bulbous lobes superimposed 
onto a series of concentric arcs, which may be associated 
with spherically-symmetric AGB mass loss (Su et el. 1998; 
2003). The lobes themselves are not uniformly illuminated, 
instead exhibiting an “S” shaped geometry (Ueta et al. 
2000). IRAS 17441-2411 is probably orientated close to the 
plane of the sky, since the object has a well-defined dark lane 
across its waist which obscures the central star. The star is 
itself extremely red; V-K>11.5 (Su et al. 2003). 

Weintraub et al. (1998) report the detection of weak H 2 
emission toward the core region in IRAS 17441-2411 (Davis 
et al. 2003b failed to detect H 2 in this source, probably 
because of the increased shot noise in their lower-spectral- 
resolution observations). Weintraub et al. postulate that the 
H 2 emission may be associated with an expanding equatorial 
torus. However, their slit was orientated east-west, rather 
than along the almost north-south bipolar lobes. 
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Figure 7. PV diagrams showing the H 2 2.122 fim emission in IRAS 17423-1755 (see Figure 2 for details). The top of each spectral image 
is orientated roughly to the north-west (p.a. ~ —45°); the “left” and “right” spectra were offset 1.0^^ to the north-east and south-west, 
respectively. The contour levels measure 0.02,0.04,0.08,0.12,0.16,0.32 Jy per arcsec along the slit in all four plots. 
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Figure 8. PV diagrams showing the H 2 2.122pm emission in IRAS 17441-2411 (see Figure 2 for details). The top of each spectral 
image is orientated roughly to the north (p.a. 17°); the “left” and “right” spectra were offset 1.2" to the east and 0.8" to the west, 

respectively. The contour levels measure 0.0025,0.005,0.0075,0.01,0.02,0.04,0.08 Jy per arcsec along the slit in all four plots. 
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Figure 9. PV diagrams showing the H 2 2.122 fim emission in IRAS 17150-3224 (see Figure 2 for details). The top of each spectral image 
is orientated roughly to the north-west (p.a. ~ —55°); the “left” and “right” spectra were offset 0.8” to the north-east and south-west, 
respectively. The contour levels measure 0.0015,0.0025,0.005,0.01,0.02,0.04,0.08,0.16,0.32 Jy per arcsec along the slit in all four plots. 


In Figure we again show continunm-subtracted H 2 
spectral images at three slit positions across the width of the 
nebula. As with the other PPN, the H 2 emission is confined 
to discrete knots in the bipolar lobes; in IRAS 17441-2411 
these are offset approximately l”-2” from the central contin- 
unm peak in our PV diagrams. In optical images the lobes 
appear to be somewhat limb-brightened. In H 2 we detect the 
strongest line emission along the central axis, although no¬ 
tably the northern, red-shifted knot (positive offsets in Fig¬ 
ure n is brighter in the “left” slit than the emission feature 
in the counter-lobe, with the inverse seen in the “right”slit. 
Note also that the velocity difference between the blue and 
red lobes is marginally larger in the left and right slits than 
it is along the central slit. Curiously, this implies higher ra¬ 
dial velocities along the edges of the bipolar lobes than along 
the central axis. Higher spatial- and spectral-resolution ob¬ 
servations would be useful to confirm this finding. 

Overall, the low radial velocities in H 2 support the idea 
that the system is orientated close to the plane of the sky. 
Emission lines in H 2 spectra extracted from the two bipo¬ 
lar lobes are again very narrow; Gaussian fits yield FWHM 
values of ^ 25 km s“^. 

3.5 IRAS 17150-3224 

IRAS 17150-3224 (AFGL 6815) also appears to be orien¬ 
tated close to the plane of the sky, with its equatorial 
torus viewed nearly edge-on (Kwok, Su & Hrivnak 1998). 
The central G-type star is largely obscured at optical wave¬ 
lengths and is highly reddened (V-K>11.0; Su et al. 2003). 
Like IRAS 17441-2411, in optical and near-IR HST images 
a series of concentric arcs are superimposed onto the bipo¬ 


lar lobes; these may again be associated with spherically- 
symmetric mass loss from the AGB phase (Kwok et al. 1998; 
Su et al. 2003). Overall, the angular extent of the PPN, 
measured in the optical across both bipolar lobes, is about 
8”- 10 ”, although fainter rings of emission extend radially 
further outward. 

H 2 emission was detected by Weintraub et al. (1998) 
and Davis et al. (2003b) in both bipolar lobes. Our PV di¬ 
agrams in Figure]^ reveal the kinematics of these emission 
features. The plots show that the emission extends across the 
width of the nebula, and that once again it is confined to 
knots in the two bipolar lobes. The knots are offset by ^2" 
from the central continuum. The emission is therefore ex¬ 
cited roughly half-way along the optical bipolar lobes, rather 
than at the ends of the lobes. 

Much lower radial velocities are recorded in 
IRAS 17150-3224 (Vlsr - 10 - 12 km s-\ FWHM 
~ 35 km s“^) than were seen in some of the other targets, 
as one would expect for a system orientated close to the 
plane of the sky. (The same peak velocity was measured by 
Weintraub et al. (1998) in their R ~ 65,000 H 2 2.122/rm 
spectrum obtained with an east-west slit.) 

3.6 H 2 non-detections in IRAS 09371-|-1212 and 
IRAS 17106-3046 

Two other PPN were also observed at 2.12/rm, 
IRAS09371-kl212 (Frosty Leo) and IRAS 17106-3046 
(see Table 1). Both targets appear to be bipolar PPN, much 
like those described above, and both are associated with 
collimated, knotty outflows (e.g. Langill, Kwok & Hrivnak 
1994; Sahai et al. 2000; Kwok, Hrivnak & Su 2000). Even 
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Table 2. H 2 2.122 fim and Br 7 line-profile fitting results from the spectra in Figure^^ Peak velocities and FWHM widths were obtained 
from either Lorentzian or Gaussian fits. The former gave a better fit to the Br 7 data; the latter more closely matched the H 2 profiles. 


Target 

H 2 (blue) 
^^eak 
(km s“^) 

H 2 (blue) 
FWHM 
(km s“^) 

H 2 (red) 
^^eak 
(km s~^) 

H 2 (red) 
FWHM 
(km s“^) 

Bry 

Ipeak 

(km S“^) 

Bry 
FWHM 
(km s“^) 

IRAS 21282-1-5050 

14.5 

27.6 

14.7 

29.7 

26.8 

25.8 

M 1-92 

-33.8 

43.3 

28.7 

37.2 

-11.4 

87.4 

IRAS 17423-1755 

-69.7 

17.7 

- 

- 

43.5 

83.7 

IRAS 17441-2411 

-1.0 

25.9 

5.9 

26.4 

4.9* 

65.3* 

IRAS 17150-3224 

10.9 

36.3 

12.3 

34.6 

20.4* 

93.1* 


* HI lines seen in absorption 


so, H 2 2.122/rm line emission was not detected along their 
central bipolar axes. 

The angular extent of IRAS 09371-1-1212 on the sky is 
considerably larger than the other targets in our sample, so 
it is certainly possible that our chosen slit position angle 
missed line-emission features along the walls of the nebula. 
Our slit was aligned through bright, optical ansae at the 
ends of the bipolar lobes (labeled An and As by Sahai et al. 
(2000)), features which tend to be H 2 line emitters in other 
PPN. Also, the outflow itself may have changed direction 
in recent years, as evidenced by the multiple, miss-aligned 
optical features observed in the lobes of this nebula (Sahai 
et al. 2000). 

IRAS 17106-3046 has a much smaller angular extent on 
the sky, though it is still comparable in size to IRAS 17423- 
1755 and IRAS 21282-1-5050, PPN where H 2 was detected 
and spatially resolved. 

Inclination angle does not seem to be a factor in whether 
H 2 is detected or not. In the other targets discussed above, 
H 2 is detected in sources with a range of inclination angles 
(see Table 1). 

Both targets certainly show evidence for collimated 
flows, and the knots and ansae in IRAS 0937H-1212 in 
particular indicate the presence of shocks. Shock velocities 
should be high enough to excite H 2 into emission; the l-OS(l) 
line targeted here can be excited in shocks with velocities of 
^ 10 km s“^ (e.g. Smith 1991). The lack of detectable emis¬ 
sion must therefore be related to a lack of dense, molecular 
gas in the polar lobes. 

Our H 2 non-detections in IRAS 09371-1-1212 and 
IRAS 17106-3046 may be related to their progenitor mass 
and/or evolutionary state. Weintraub et al (1998), Garcfa- 
Hernandez et al. (2002) and Davis et al. (2003b) all find that 
H 2 is not detected in targets with spectral types later than 
about G2. Our non-detection in IRAS 09371-1-1212, which 
has a K7II central object (Mauron, Le Borgne & Picquette 
1989), would support this hypothesis. The spectral type of 
IRAS 17106-3046 is not known. 


3.7 Bry observations 

Bry observations were acquired for the central slit position, 
using the same slit position angle, in each of the five targets 
where H 2 was detected. In two of the five sources - the two 
late, G-type sources in our sample - Bry is seen in absorp¬ 
tion. This is typically the case in F or G-type PPN (e.g. 


IRAS21282+5050 



Velocity (km/s) 
MI-92 



Velocity (km/s) 
IRAS1 7423-1755 



Figure 10. PV diagrams showing the Bry emission in three of 
the five targets observed (Bry was seen only in absorption in 
IRAS 17150-3224 and IRAS 17441-2411). The continuum emission 
(at zero offset in the y-axis) has been fitted and removed. In 
all three plots the contour levels start at 0.025 Jy arcsec“^ and 
increase in multiples of 2. 

Hrivnak, Kwok & Geballe 1994). In the “more evolved” B- 
type PPN and the “young PN” IRAS 21282-1-5050, Bry is 
strong and seen in emission, indicating the onset of ioni¬ 
sation. Ml-92, IRAS 17423-1755 and IRAS21282-h5050 are 
probably on the brink of the transition from PPN to PN 
(Kwok 1993). In these three targets the emission is largely 
confined to the source continuum position fFigure ITDll . al¬ 
though weak, blue-shifted emission, presumably associated 
with the inner regions of the post-AGB wind, is also seen in 
each case. 

In Figure EU we compare Bry profiles to H 2 2.122 /rm 
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Figure 11. H 2 2.122 /im and Br 7 spectra: H 2 spectra have been 
extracted over a 2.7"-wide strip ~ 2" above and below the source 
continuum and are drawn with dashed lines (in IRAS 17423-1755 
H 2 emission was only detected in one bipolar lobe); the Bry spec¬ 
tra are extracted over the same area though these are centred on 
the source continuum position. The H 2 spectra for Ml-92 have 
again been scaled by a factor of 0.1. 


profiles, the latter being extracted from the bipolar lobes of 
each PPN. The central slit images were used for the H 2 spec¬ 
tra. The results of Gaussian (H 2 ) or Lorentzian (Br 7 ) fits to 
these lines are listed in Table 2. The Br 7 profiles generally 
peak at low radial velocities, velocities that are compara¬ 
ble (to within ±10 km s“^) to systemic velocities derived 
from CO observations of the molecular gas associated with 
each system (Bujarrabal et al. 20001 and references therein). 
Considering the targets individually; 

• The two H 2 profiles from the young O-type PN 
IRAS 21282±5050 peak at essentially the same low veloc¬ 
ity, as one would expect for a system orientated close to the 
plane of the sky. The Bry profile peaks at a similarly low 
velocity, and is notably much narrower than the Bry profiles 
observed toward the other sources (Table 2). 

• In the B-type PPN Ml-92, the Bry emission peaks mid¬ 
way between the two H 2 emission lines. Here, the HI recom¬ 
bination line is probably centred at the systemic velocity, 
while the H 2 lines reflect the blue- and red-shifted gas ve¬ 
locities in the bipolar lobes. 

• In the B-type PPN IRAS 17423-1755 there is a consid¬ 
erable shift between the H 2 and Bry peak velocities. Again, 
the Bry line peaks at the systemic velocity (as reported by 
Riera et al. 2003), while the H 2 traces the high-velocity blue- 
shifted outflow lobe. The red-shifted lobe may be obscured 
by a circumstellar torus inclined toward the plane of the 
sky. The large velocity shift 130 km s“^) between the 
H 2 and Bry profiles would likewise be caused by the large 
inclination angle of the system. 

• In the G-type PPN IRAS 17441-2411 and IRAS 17150- 
3224 the Bry absorption features again probably mark the 
systemic LSR velocity in each system. The H 2 profiles in 
the bipolar lobes in both targets peak at roughly the same 
radial velocity, as one would expect for nebulae orientated 
close to the plane of the sky. 

In four of the five sources (the PPN) the extracted Bry 
profiles are considerably broader than the H 2 profiles, by 
a factor of ^ 2 (Table 2). The Bry emission must trace 
ionised regions within the fast post-AGB wind, while the 
H 2 instead reflects the shock velocity at the interface be¬ 
tween the ionised wind and the slow-moving AGB envelope. 
Note in particular the broad blue and red-shifted wings su¬ 
perimposed onto the otherwise Lorentzian profile in M 1-92. 
These wings reach velocities of ±200 km s“^ (at the 10% 
peak intensity level). For an inclination angle of ~ 35° this 
indicates a post-AGB ionised wind speed of at least 350 km 
s“^. 

In contrast with the other sources, the Bry profile to¬ 
ward the young PN IRAS 21282±5050 is very narrow, being 
comparable in width to the H 2 profiles. We also see no ev¬ 
idence for associated line wings. Together, the H 2 and Bry 
data support Meixner et al (1997)’s claim that the bipolar 
lobes are orientated close to the plane. Our data also suggest 
that the HI recombination line emission must be excited in 
a poloidal rather than equatorial wind, and that this wind 
must be highly collimated, perhaps by the torus reported 
by Meixner et ah, since any lateral expansion would serve 
to broaden the line profile. 









12 C.J. Davis et al. 


4 FURTHER DISCUSSION 

4.1 Br 7 and H 2 as tracers of the fast wind 

Br 7 and H 2 clearly trace different components in the out¬ 
flows from the sources considered in this paper. Br 7 emission 
is produced close to the central continuum source, the line 
profiles are generally broad though they peak close to the 
systemic velocity, and the Bry flux in each system rapidly 
decreases with distance along the bipolar axes. H 2 , on the 
other hand, is predominantly excited in the bipolar lobes 
in each PPN and appears to be confined to shocks within 
each flow (we detect molecular line emission toward the 
continuum-peak in only one system, IRAS 21282-1-5050). H 2 
profiles are narrow, though their peak velocities are blue- 
and red-shifted in bipolar nebulae that are inclined out of the 
plane of the sky (note in particular IRAS 17423-1755). We 
will discuss the origin of the H 2 emission in subsequent sec¬ 
tions. However, we first consider possible excitation mecha¬ 
nisms for the observed Bry emission. 

Detectable Br 7 recombination-line emission requires 
high gas densities or an intense ionising flux. The former 
is produced in post-shock gas; the latter exists close to in¬ 
termediate or high-mass stars. High-resolution imaging and 
kinematic studies of PPN show them to be associated with 
outflows and shock phenomena. PPN, as precursors to PN, 
are also potentially associated with compact, photo-ionised 
regions. Either environment may therefore be a source of the 
observed Brq emission. 

Atomic hydrogen is collisionally-ionised in shocks, pro¬ 
vided the post-shock gas temperature exceeds ^ 4 x 10'*K, 
i.e. for shock velocities greater than about 50 km s“^. Radia¬ 
tive recombination then leads to bright Lyman and Balmer 
emission. The higher-energy Bracket lines will be much 
weaker, of the order of 0.01-0.03 times the H/3 intensity for 
shock velocities up to 150 km s”*^ and gas densities in the 
range 10® — 10^cm“® (Hollenbach & McKee 1989). Indeed, 
for a shock velocity of ^100-150 km s“^ and a pre-shock gas 
density of ~ 10"^ cm“®, Hollenbach & McKee predict a Bry 
line intensity of the order of 10“®'^ erg cm“^ s“® sr“^. For 
a shock covering an area of 1” (2.35 x 10“^^ sr) on the sky, 
and for a line width of ~ 80 km s“^ (Table 2), this is equiv¬ 
alent to a Bry line flux density of about 0.01 Jy. This is an 
order of magnitude lower than the observed flux levels in the 
two early-type PPN and the young PN, IRAS 21282-1-5050 
IFigure im . Much higher gas densities, of the order of 10®- 
10® cm“®, would be needed to produce the observed line 
strengths. 

In an ionised region at the base of a fast wind, Bry 
will again be relatively weak in comparison to Lyman and 
Balmer lines, ^0.03 times the H/3 intensity for an electron 
temperature and density of ~ 10“^ K and ~ 10"^ cm“® (un¬ 
der Case B [optically thick to Lya] recombination theory; 
Hummer & Storey 1987). The total integrated intensity in 
Bry from an optically thin PPN nebula of radius i? at a 
distance D will be (Kwok 2000): 

/sr-7 = 3.41 X lO’^’^neTip ^rg cm“®s“^. 

For D = 2kpc, R = 0.01 pc (~2”), rie ~ rip = 10"^ cm“® 
and a filling factor, e, of unity, the total Bry flux will be 
~ 0.9 X 10“^® erg cm“^ s“^. For a line width of ~ 80km 


s“® (Table 2), the flux density will be ~ 0.2 Jy, comparable 
to the values in our extracted spectra IFigure lTTll . 

Thus, from the observed line intensities, a photo-ionised 
nebula close to the central object seems to be the most likely 
origin for the Bry emission. The spatial distribution of Bry 
along the slit also tends to support recombination at the 
base of the outflow, rather than in shocks between the fast 
wind and the AGB envelope, since these shocks, traced in 
H 2 in each system, are at a resolvable distance from the 
central source and yet the Bry emission is largely confined 
to regions near the central source. The low peak velocities 
though broad line widths in the Bry spectra also support 
recombination in a fast, possible poorly collimated or even 
spherical wind close to the source. 

Finally, we mention that the observations of Guerrero et 
al. (2000) and Davis et al (2003b) and the models of Natta & 
Hollenbach (1998) indicate that the H 2 2.122 /rm/Bry ratio 
is highest in PN with large, well-defined rings and bright cen¬ 
tral stars, i.e. in evolved and/or intermediate-to-high mass 
sources. The ratio seems to be lower in “less-evolved” PN. 
When we consider PPN, particularly those with late-type 
central stars like IRAS 17441-2411 and IRAS 17150-3224, 
where thermal and probably shock-excitation dominates the 
production of H 2 line emission, the H 2 2.122/im/Bry ratio 
is again high. 

4.2 The “three wind” model 

Three phases of mass loss have been invoked to explain the 
multiple shells observed in some PN (e.g. Frank 1994). In 
this scenario the phase of slow, spherical AGB mass loss 
(M~ 10“^-10“® Mq yr“®; U ~ 10 km s“®) is followed by a 
slow though more massive “superwind” (M~ 10“®-10“® Mq 
yr“^; U ~ 10 km s“^). A diffuse, fast wind (M~ 10“® Mq 
yr“^; V ~ 1000 km s“^) follows a few thousand years after 
the superwind has been terminated. The model is probably a 
simplification of a continually changing post-AGB mass-loss 
rate, and in Frank (1994)’s treatment only the “preionisa¬ 
tion stage” is relevant to PPN. Nevertheless, we investigate 
briefly whether our near-IR spectra can be used to distin¬ 
guish between the traditional “Two-Wind” model of Kwok 
(1993) and a more complex scenario. Is an additional super¬ 
wind required to explain the H 2 observations? 

The density discontinuity between a slow superwind 
and an equally-slow AGB envelope drives a boundary shock 
(Frank 1994) which could produce strong H 2 line emission. 
However, it would not generate the high radial velocities seen 
in our H 2 spectra (velocities approaching 50-100 km s“^ in 
some case; Table 2) since neither the superwind nor the AGB 
wind expand at high speed. So the observed H 2 emission 
must not derive from a Superwind/AGB wind boundary. 
Instead, a fast wind impacting the back of either a super¬ 
wind or an AGB wind is needed to produce the observed 
high velocities. 

The boundary layer between a diffuse, fast wind and 
a dense, neutral superwind or AGB wind would comprise 
a weak, outer/forward shock and a strong, inner/reverse 
shock. The latter, with Uhock ~ 1000 km s“^, decelerates 
and heats the fast wind material to very high temperatures 
{T ~ lO^K). The velocity of the forward shock would be 
much lower, comparable to though somewhat higher than 
the velocity of the AGB wind, depending on the ratio of 
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Figure 12. A bow shock model for the H 2 emission from IRAS 21282+5050. The H 2 profiles in the third panel can be compared to the 
extracted spectra in Figure U] In the model, a parabolic bow of speed 70 km s ^ advances into a uniform medium of hydrogen density 
njl = 10^ cm“® (of which 60% is in atomic form) and Alfven speed 1 km s“^ (B ~ 50 pG), with the magnetic field parallel to the 
bow motion. The size (semi-latus rectum) of the bow is 2 X 10^® cm (~0.5” at 3 kpc). The line profiles in the third panel correspond to 
the space-integrated emission from the two slits shown in the first panel. The PV diagram in the middle panel has been calculated by 
integration through the entire transverse direction. The line profiles were convolved with a Gaussian of velocity dispersion 7.2 km s“^ 
and the PV diagram with 15 km s“^. In the left panel the star would be to the left; in the middle panel the star would be to the top. 


the fast wind density and the AGB or superwind density. 
For example, a fast wind velocity (t4ast) of 1000 km s“^ 
and a density ratio between the fast wind and AGB wind 
(pAGB/pfast) of 100 would result in an outer shock velocity, 
given by 

Khock ~ Vfast/[1 + (PAGB/Pfast)^^^], 

of ~ 90 km s“^. By comparison, a density ratio of 10,000 
between the same fast wind and a much denser superwind 
would yield a much lower shock velocity, of the order of 
10 km s“^. In both cases, because the shocked H 2 would 
more-or-less move with the shock front, the peak H 2 velocity 
in observed spectra (and the proper motions of H 2 emission 
knots), when corrected for the flow inclination angle, should 
be comparable to these shock velocities. Since moderately 
high H 2 velocities have been observed, our observations tend 
to support a lower density ratio and so, at least for the H 2 
emission regions considered here, a third, dense superwind 
component may not be required. 


4.3 H 2 excitation in bow shocks 

The morphologies and excitation of molecular hydrogen fea¬ 
tures in PPN, when combined with the kinematic data pre¬ 
sented here and elsewhere (e.g. Weinbtraub et al. 1998), 
strongly support excitation in curved or conical shocks. 
Models of similar “bow” shocks have been invoked to ex¬ 
plain line profiles and excitation maps in HH jets and molec¬ 
ular outflows from young stars, most recently the molecu¬ 
lar bow shocks observed in embedded systems (e.g. Smith 
1991, 1994; Davis et al. 1999; Smith, Khanzadyan & Davis 
2003; O’Gonnell et al. 2004). In these papers the effects of 
varying orientation angle, bow shock speed, bow shape, pre¬ 
shock density, and the Alfven speed and magnetic field ori¬ 
entation are investigated and discussed in detail. Here we 


present only the models which best fit the kinematic data 
obtained for two systems: IRAS 21282-1-5050, where complex 
and double-peaked profiles are observed, and M 1-92, where 
H 2 and [Fell] emission is detected. In both cases we find that 
a J-type bow shock fits the data reasonably well; in other 
words, C-shocks and therefore strongly magnetised winds 
are not required. 

Wherever possible we use parameters derived from 
other observations to constrain the models. 


4 . 3.1 IRAS 21282+5050 

The bipolar lobes of the IRAS 21282-1-5050 nebula appear 
to be orientated close to the plane of the sky. In H 2 we 
observe emission along the flow axis, but also toward the 
central continuum peak, where double-peaked profiles with 
a peak-to-peak separation of 22 km s“^ are observed (Fig¬ 
ure In Section 3.1. we discuss whether the emission to¬ 
ward the continuum could be excited in an equatorial wind. 
We have attempted to model such an equatorial outflow with 
bow shocks propagating along the line of sight. This might 
be considered consistent with the details displayed in Fig¬ 
ure]^ and discussed in Section 3.1; specifically, the double- 
peaked 1-0 S(l) profile plus single peaked profiles from the 
left and right-hand spectra. However, after detailed mod¬ 
eling we were unable to reproduce the peak shifts. This is 
because bow shocks tend to deflect material away from the 
shock surface rather than plough it up. Even with a slow 
non-dissociative bow shock, most of the emission arises from 
the wings where the emission peak is shifted by only a few 
km s“^. 

Instead, we more easily model the IRAS 21282-1-5050 
H 2 emission if it arises from bow shocks moving in the plane 
of the sky, as displayed in the left panel in Figure WI\ In 
this case, the deflection from the bow flanks is found to 
be sufficient to generate a twin-peaked H 2 profile from the 
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Figure 13. A bow shock model for the H 2 and [Fell] emission from M 1-92. The [Fell] PV diagram and H 2 and [Fell] line profiles can 
be compared to the observed red-shifted components shown in Figures U] andIn the model, a parabolic bow of speed 200 km s“^ 
advances into a uniform medium with hydrogen density njj = 10^ cm“® (60% in atomic form) and Alfven speed 1 km s“^ (B ~ 50/rG), 
where the field is parallel to the bow motion. The size (semi-latus rectum) of the bow is 2 X 10^® cm (>^0.6” at 2.5 kpc). In the centre 
panel, the [Fell] PV diagram has been calculated by integration through the entire transverse direction. In the right-hand panel, the 
[Fell] (full line) and H 2 (dotted line) profiles correspond to the space-integrated emission from horizontal slits passing through the bow 
apex and the bow wing, respectively. The PV diagram and line profiles have been convolved with a Gaussian of velocity dispersion 15 km 
s“^ and 12 km s“^, respectively. In the left panel the star would be to the left; in the middle panel the star would be to the top. 


central slit, with a peak-to-peak separation consistent with 
the observations in Figure 01 a single-peaked profile is also 
produced in the bow wings, corresponding to the left or right 
slit in Figure 01 The overall line widths are also roughly 
consistent with the observations, and the model profiles peak 
at the systemic velocity, as is more-or-less observed to be 
the case. The profiles have been convolved with a Gaussian 
profile to represent turbulent motions which will, in any case, 
be generated wherever there are curved shock fronts. 

In Figure lT^ we also reproduce synthetic H 2 (and [Fell]) 
images, for comparison with future high-resolution imag¬ 
ing, and an integrated PV diagram, for comparison with the 
plots in Figurel^ (Note that the model only yields emission 
from one bow shock and therefore only one outflow lobe; 
in the observed PV diagrams the emission is presumably 
excited in both lobes.) 

The model employed solves the 1-D MHD equations 
at each point on the bow surface. The cooling behind the 
shock and the motion away from the bow surface is taken 
into account with the assumption that the cooling length 
is short (but not negligible) in comparison to the bow size. 
The orientation derived from this model is consistent with 
the direction of the bipolar outflow on larger scales (Table 1). 
However, we cannot exclude other models given the present 
limited dataset and modest resolution. 


4 . 3.8 Ml-92 

M 1-92 was the only PPN where we detected H 2 and [Fell] 
emission: we therefore attempted to model both emission 
lines, using bow shocks orientated at various angles to the 
line of sight. We also varied the density, bow shape, field 
strength, molecular fraction and bow speed. We obtained 
good agreement with the data for a high-velocity bow shock 
propagating at 35±10° to the sky plane, for a density ex¬ 


ceeding 10® cm“® and a bow shape s ~ 2.0 ± 0.2, where the 
bow shape is defined hj R oc Z’’ in cylindrical coordinates 
(e.g. Smith et al. 2003). The density, velocity and inclina¬ 
tion angle used are consistent with those derived through 
modeling of atomic lines in the optical (Bujarrabal et al. 
1998a). 

The model results are shown in Figure^! where we 
again display synthetic H 2 and [Fell] images of the bow 
structure. As is typically observed in HH bow shocks, the 
[Fell] emission is confined to the apex of the bow, with the 
H 2 brightest in the non-dissociating bow flanks. Because the 
M 1-92 bow model is advancing at a much higher speed than 
the IRAS 21282-1-5050 bow, the [Fell] emission is more ex¬ 
tended and the H 2 is pushed further back into the bow flanks 
(compare Figure^] with Fieure lT^ . In high-resolution H 2 
images of M 1-92, the H 2 emission does seem to be confined 
to off-axis peaks that could be the wings of a bow shock 
(Bujarrabal et al. 1998a; Davis et al. 2003b). From profiles 
extracted along the slit axis (Section 3.2), we also find that 
the [Fell] peaks in the two lobes of the bipolar nebula are 
slightly further apart than their H 2 counterparts, as would 
be expect from two oppositely-directed bow shocks. 

The predicted PV diagram for the [Fell] data is shown 
in the second panel of Figure ESI The plot displays a boot¬ 
shaped structure comprising an extended red wing and faint, 
low-velocity emission trailing back toward the source. The 
overall shape of the model PV diagram agrees quite well with 
the observed profile in Figure El for the southern emission 
lobe, although the model does predict double-peaked profiles 
which are not observed. 

Extracted H 2 (dotted line) and [Fell] (full line) spectra 
are shown in the third panel in Figure El The model repro¬ 
duces the general structure of both the H 2 and [Fe II] lines; 
i.e. the narrow, Gaussian H 2 profile from the bow wings to¬ 
gether with the broader, asymmetric [Fell] profile from the 
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bow apex (compare the dotted H 2 and full-line [Fell] model 
spectra in Figure IT!^ with the observed spectra in Figure]^. 

The extreme [Fell] line wings, which extend out to > 
300 km s“^, are not reproduced by the model. Moreover, 
the predicted H 2 and [Fell] lines peak at or very near to 
the velocity of the ambient medium, which was taken to be 
stationary. The observed profiles are red-shifted (and blue- 
shifted) by 35 — 45 km s“^ (Figure]^. The observations 
imply that the pre-shock medium is expanding with a radial 
speed of ~70 km s“^. Such a fast-moving pre-shock medium 
ahead of each bow would be sufficient to shift the model line 
profiles and PV diagram by about 35 km s“^. 

As with IRAS 21282-1-5050, the bow models do an ade¬ 
quate job of explaining the observations, although additional 
high-spatial-resolution near-IR imaging and spectroscsopy 
would again be useful to more tightly constrain the theory. 

4.4 Rotation in a magneto-centrifugal disk wind 

Finally, since we have observed H 2 spectra at three positions 
across the width of each PPN, we are potentially in a posi¬ 
tion to search for signs of rotation. Hints of jet rotation have 
recently been identified on large scales in an outflow from a 
young star by Davis et al. (2000), and investigated in more 
detail in HH jets by Bacciotti et al. (2002) and Coffey et 
al. (2004). In magneto-centrifugal disk winds, the outflow is 
expected to maintain a record of disk rotation, particularly 
close to the central star, so their findings are not unexpected. 

IRAS 17150-3224 is the only source that shows possi¬ 
ble signs of rotation. This source is orientated close to the 
plane of the sky, is one of the brightest H 2 emitters in our 
sample, and has emission features well resolved from the 
central continuum. In Figure El we see a trend of increasing 
velocity from right to left across the width of the outflow. In 
other words, the south-western edges of both lobes appear 
to be moving toward the observer, with the north-eastern 
edges moving away from the observer, with a toroidal veloc¬ 
ity of 3-7 km s“^, depending on which lobe is considered. 
The trend is consistent across all three slits, and is in the 
same sense in both lobes. If we adopt the slit offset (0.8", 
or 3 X m at a distance of 2.4 kpc) as the radial dis¬ 
tance, then an angular velocity of ~ 2 x rad s“^ and 

an orbital period of 20,000 years is implied. 

Disk wind models are reviewed by Konigl & Pudritz 
(2002) and discussed in the context of jet rotation by Bac¬ 
ciotti et al. (2000) and in terms of PPN by Frank & Black¬ 
man (2004). Briefly, in a disk wind material is ejected from 
the surface of the disk at a “footprint” radius, ro, and carried 
poloidally outward along field lines, provided the base of the 
field lines are inclined at less that 60° to the disk plane. The 
point along the flow where the poloidal flow velocity reaches 
the poloidal Alfven speed is known as the Alfven surface; 
the flow radius at this point is called the Alfven radius. The 
outflowing gas is accelerated only between the disk surface 
and the Alfven surface, where gas motions can be approx¬ 
imated to solid body rotation. In this region the toroidal 
velocity is proportional to the radius; v^f, oc r. Beyond the 
Alfven surface angular momentum is conserved as the gas 
moves along the field lines; in this region v^j, oc r~^. Conse¬ 
quently, at some point further downwind along the flow axis, 
the observed linear rotational velocity and radius, Uoba and 
robs, can be related to the rotational velocity and radius at 


the Alfven surface, va and ta, and the Keplerian velocity 
and radius at the flow footprint, Vo and r^, by: 

Wobsrobs ~ VArA ~ Uo(rA/ro)rA. (1) 

We assume that the flow does not interact with (and be ap¬ 
preciably decelerated by) the ambient medium: hence, mea¬ 
surements as close as possible to the disk are desirable. For a 
“cold” flow the maximum poloidal velocity along the mag¬ 
netic field lines will be equivalent to the outflow velocity. 
This is given by: 

Wflow ~ 2^/^Uo(rA/ro). (2) 

With measurements of Uobs and robs, together with an esti¬ 
mate of the outflow velocity uhow, from equations 1 and 2 
one can estimate the Alfven radius, ta, i.e.: 

ta ^ 2 ^ Uobsrobs/rflow. (3) 

Disk wind models typically predict a value of 2-5 for the 
ratio rA/ro. If we assume a value of 3, then an estimate of 
rA yields an estimate for the footprint radius, ro. 

Finally, from Kepler’s law (uo = {GM/roY^'^), we may 
estimate the mass of the central star: 

GM ~ 2“^^^Wobsrobsrflow(ro/rA)^. (4) 

In this way the footprint radius for the disk wind and the 

mass of the central object may be derived using observable 
qualities, Uobs, robs and uhow, provided the theory is correct. 
For IRAS 17150-3224 our slit spectra imply a toroidal ve¬ 
locity Hobs ~ 5 km s“^ at a radius robs ^ 3 x m. The 
poloidal velocity of the molecular flow traced in H 2 is diffi¬ 
cult to measure because the bipolar nebula lies close to the 
plane of the sky. We assume a canonical speed of 100 km s“^ 
which, judging from our observations of other PPN, does not 
seem unreasonable. These values lead to an Alfven radius of 
^ 134 AU and a source mass of 30 M©. 

Both seem to be somewhat on the high side (e.g. Frank 
& Blackmann 2004; Kwok 1993). Because we are observing 
the molecular component of the flow, a lower flow velocity 
might be considered more approriate. This would serve to 
reduce M, though va (and ro) would then increase. Large 
Alfven and footprint radii might be expected for the molec¬ 
ular component of a disk wind, since molecular hydrogen 
would probably not survive acceleration close to the central 
object. Alternatively, robs could be overestimated, particu¬ 
larly given the modest spatial resolution of our observations, 
and the fact that the bow shock will be broader than the 
underlying outflow. M and va are both directly proportional 
to the observed flow radius; reducing robs would lower both 
quantities. 

The interaction of the fast wind with the AGB envelope 
could of course produce asymmetries across the width of the 
flow. In particular, entrainment along one side of the flow 
could decelerate gas along the flow edge producing skewed 
velocities that mimic rotation. Such an effect would not be 
expected in both lobes of the flow, however, and proba¬ 
bly would not explain the gradual increase in velocity seen 
across the three slit positions. 

We conclude that, give the limited spatial resolution 
afforded by our data, IRAS 17150-3224 is certainly a candi¬ 
date for nebula rotation, one that deserves further study in 
a number of emission lines with a high-resolution adaptive- 
optics spectrometer. 
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5 CONCLUSIONS 

From our high spectral-resolution near-IR spectroscopy we 
arrive at the following conclusions: 

(i) In bipolar PPN, H 2 emission is excited in - and largely 
confined to - molecular shocks in the bipolar lobes. The 
emission is probably excited in shocks between the slow- 
moving AGB ejecta and the much faster, post-AGB wind, 
specifically in the dense, swept-up layer behind the forward 
shock at the fast wind/AGB wind interface (Frank 1994). 
H 2 velocities are certainly more comparable to the velocity 
of the detached AGB envelope than the fast wind. 

(ii) Unlike Garci'a-Hernandez et al (2002), we find (albeit 
in a very limited sample biased toward bipolar PPN) bright 
H 2 emission in PPN associated with early and reasonably 
late-type central stars. Our sample includes one O, two B 
and two G-type objects. However, we did fail to detect H 2 
along the axis of a bipolar PPN associated with an M-type 
star. 

(iii) H 2 emission is associated with PPN, regardless of 
whether they exhibit Bry in emission or absorption. This is, 
at least in part, due to the fact that the H 2 and HI regions 
are physically separate. 

(iv) Our H 2 observations of IRAS 21282-1-5050, and both 
H 2 and [Fell] observations of M 1-92, are convincingly mod¬ 
eled with molecular bow shocks viewed roughly side-on. Line 
widths, peak velocities, and double-peaked profiles are re¬ 
produced. The models also predict a spatial offset between 
the H 2 and [Fell] emission regions, the former being excited 
in the oblique bow flanks, the latter at the bow apex. There 
is evidence for such an offset in M 1-92, although additional 
observations at higher spatial resolution in these and other 
PPN are warranted. 

(v) We find possible evidence for outflow rotation in one 
target, IRAS 17150-3224. Rotation is expected if the flow 
is driven magneto-centrifugally from the surface of a disk. 
Spectroscopy at high spatial and spectral resolution would 
be extremely useful, with a system that operates in the near- 
IR (such as NACO at the VLT) to confirm our results and 
to search for similar trends in other targets, since rotation 
would validate disk wind models, which have recently been 
invoked to explain other aspects of outflows from evolved 
stars. 
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